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We present a self consistent treatment of pinhole junctions in d x 2_ y 2 superconductors. The 
current-phase relation j a (x) is studied at different temperatures and at different angles a between 
the crystal a-axis and the junction (surface) normal. We show that the critical current of a junction 
can be reduced by pair-breaking effects at the separation. We also study the Josephson energy 
of a pinhole as a function of the phase difference across the junction. In particular are mapped 
the positions of the energy minima x-min at different temperatures as functions of {ctL,a>R), the 
crystal orientations of the left and right superconductors. With decreasing temperature there is an 
increasing range of crystal orientations where Xmin varies continuously from to n. 
PACS numbers: 74.50.+r, 74.72.-h 
(— i | keywords: superconducting junctions, unconventional pairing 
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Recent experiments and their interpretations have strongk argued in favor of an order parameter of B\ g ( d x 2_ y 2 
) symmetry in-the oxide superconductors (see e.g. reviewsua). Part of the evidence comes from bulk and transport 
properties (seeEl). Other relevant data involve Josephson junctions which probe the order parameter at surfaces. In 
these latter experiments, the phase of the order parameter is involved, and "interference" experiments can be carried 
out. 

It will be important below to realize that interference experiments fall into two categories. In the first one measures 
the critical current of a two junction loop as a function of the magnetic flux treading the ring and records the 
positions of the current maximaflu. The second type consists of mapping the phase differences at the Josephson 
energy minimaLTLj. 

Josephson junctions in s-wave superconductors have been well explopd. Here we distinguish between "tunnel" 
junctions on the one hand and "weak-links" on the other. We shall followH-3 and classify junctions on the basis of their 
normal-state conductivities. In a tunnel junction normal electrons face a reasonably high potential barrier whereas in 

■ a weak-link junction they do not. 

y—i | Even in s-wave superconductors, tunnel junctions and weak links behave in different fashions. We compare weak- 
Q\ . link junctions with tunnel junctions in unconventional superconductors picking the pinhole as the representative of the 
former. The differences arise from two sources: Firstly, the current-phase relation of a tunnel junction is sinusoidal. 
The minima/maxima of the Josephson energy therefore lie at a phase difference equal to zero or n (provided the order 
paramete r do es not break time-reversal symmetryllj) . In a weak link, the current-phase relationship is not sinusoidal 
in generaOO, and the minima/maxima will not lie at zero or it. Secondly, in weak links the current carrying processes 
sample a wider angular range of momentum than in tunnel junctions. With different crystal orientations on the two 
sides of a junction, varying the angle of the junction with respect to the crystals, one observes quite different behaviors 

■ of the positions of the energy minima and critical current in pinholes and in tunnel junctions. 

Nevertheless, " weak- links" , (especially our pinhole) and tunnel junctions are the opposite ends of a continuum. 
Imagining a weak potential barrier in the link and increasing its strength should take one smoothly from the pinhole 
' to the tunnel junction. 

^ , JosepJaspa junctions are probes of surfaces. Near a surface, an unconventional order parameter must be strongly 
affectecOH We therefore also investigate the effect of surface pair breaking on our weak links and tunnel junctions. 

Most of the discussions on high-Tc junctions in the literature ignore the difference between weak links and tunnel 
junctions (and the effects of surface pair- breaking) . In the sequel it is argued in some detail that a number of 
puzzling reports in the literature may have a natural explanation if one keeps track of these effects. One should 
be very careful when drawing conclusions on the order parameter type of a high T c superconductor on the basis of 
measurements of various Josephson-related effects. All available models and their variants should checked out for 
competing interpretations. 
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In Section 2 we discuss the self-consistently computed junctions, first pinholes (weak links) and then tunnel junctions 
in Section 3. Section 4 is about the non self-consistent junctions, i.e. those treated analytically assuming constant 
order parameters up to the separation where the junction sits. Section 4 ends with the Subsection 4.1 which compares 
the self-consistent and non self-consistent models introduced in the previous Sections. Section 5 is on critical currents 
and phase differences where the Josephson energy reaches its minima. The final Section 6 is on the conclusions we 
suggest from experiments. The ultimate Section 6.3 is a critical analysis of the plausibility of our models as describing 
real junctions. [] 

II. THE PINHOLE JUNCTION 

The pinhole junction is a small opening in an interface separating two superconductors (see Figure 0), small both 
in length and in width on the scale of the coherence length £o- Being so tiny, the opening perturbs little the order 
parameter which may be calculated ignoring the pinhole. A phase difference x over the junction is introduced 
multiplying the self-consistent order parameters on the two sides with the factors exp(±zx/2). The phase is,taken to 
hop discontinuously at the partition. The pinhole junction was first introduced in s-wave superconductorsc2l and has 
also been considered as a model of a weak link in superfluid 3 HeEilt3. 
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FIG. 1. Sketch of a pinhole. The opening is small on the scale of coherence length. Most trajectories are reflected at the 
interface, and only those hitting the hole are transmitted. 

A. The self-consistent pinhole 

Surfaces have a pair-breaking effect on a superconductor whose order parameter depends on the direction p of 
propagation of a quasiparticle, and we need to take into-,account the quasiparticle scattering at the walL-ModelsJor 
walls have been studied in 3 He (see, for instance^ 3 ! andEa) and recently also in d-pWjave superconductors (Il3j2a andc3). 
We chose the specular surface, attractive in its simplicity. It has been studied bjfia. The Green's function g(p, R; e„) 
in the specular model is taken to be continuous at the surface along pairs of trajectories (pin, Pout) connected as 

Pout = Pin - 2n(n • p m ), (1) 

n being the surface normal. 

With the boundary conditions set, the procedure of getting at the order parameter consists of iterating the Eilen- 
berger equation 



x We shall only consider planar interfaces. Some recent papers ( e.g££l and references therein) have suggested that it may be 
important to include the effects of facets in understanding the properties of some grain boundaries. If this is so, then we cannot 
directly apply our results to these grain boundaries; instead they can only serve as inputs when considering the facets. 
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[ie n f 3 - A(p, R), g(p, R; e n )\ + iv f ■ V R <?(p, R; e„) = (2) 

for the propagator, where Vf is the Fermi velocity and p the direction of momentum at the Fermi surface, together 
with the gap equation for the order parameter 

A(p,R)=7rT£ J ^lF( p ,p')iTr{(n-if 2 )5(p',R;e n )} (3) 

till self consistency. As indicated by the carets, the propagator g and the order parameter A are matrices in the 
particle-hole space The order parameter is conveniently parameterized in terms of a real part A2 and an imaginary 
part Ai as A = i(Ai fx + A2 fa). Only interested in an equilibrium property, we can use the Matsubara formalism at 
the Matsubara frequencies e n — %T(2n + 1). 

The pairing interaction V(p,p') in equation (|J) determines the symmetry of the order parameter. V(p, p') can be 
separated into a sum of allowed pairing channels. The strength of a particular channel X depends on the coupling 
parameter V ^ X ) which can be eliminated in favor of the transition temperature T C (X) in the channel with the aid of 
the well-known BCS-relation 

—, — r = hi — ; — - + r (4) 

V(X) T C (X) $£n+\ { > 

where n c is a cutoff. In thisjsfork, the dominant pairing channel is always chosen as having an order parameter of B\ g 
symmetry. Buchholtz et alt3 have shown that for equal VlX) in the subdominant channels, i?2g has the strongest 
effect on the properties of a smooth surface. Thus we followta and study a superconductor with a d x 2_ y 2 or B\ g order 
parameter in the bulk and allow an admixture of a d xy or £?2 ff component close to the surface. The degree of mixing 
is controlled by the ratio of the transition temperatures of the two representations. This ratio is a parameter of the 
present calculation. The pairing interaction is 

V(p,p') = 2V Blg cos 2(f> cos 2cj>' + 2V B2g sin 20 sin 20'. (5) 

We refer tolll for a full account of the behavior of the order parameter at a surface and only briefly summarize the 
principal effects here. Imagine a superconducting crystal with a d x 2_ y 2 order parameter cut with a surface whose 
normal is n. In a thin film sample, n would lie normal to the junction in the plane of the film. The d-wave order 
parameter is reduced when the crystal is rotated relative the surface (junction) normal n to a position determined 
by the angle a between the crystal a-axis and the normal n (see Fig. |J The order parameter is unaffected by the 
presence of the surface (=cut) if the crystal has its a or 6-axis parallel to n. It is maximally reduced when a is equal 
to j + n^, n being an integer. In the latter case, the order parameter vanishes identically at the wall. If the ratio 
of transition temperatures is chosen larger than zero, a Bi y component may develop in the vicinity of the wall. A 
special case is at a close to \ at temperatures T /T c (Bi g ) < T c (B-2 g ) /T c (Bi g ) where the composite order parameter 
is found to break time-reversal symmetrycJ'cH. 

There is, however, a problem associated with mixing two representations as described above. On the basis of the 
boundary condition g(p>i n , 0; e„) = g(p ut, 0; e„), eq (J3J) and the general symmetries of the propagator^ one can show 
that (see alsonj), for T c (£?2 g ) different from zero, the ratio of the two order parameter components at the wall is^| 

AB - (0) tan 2c, (6) 



A Bl9 (0) V Blg 



Numerically this becomes 



2 The propagator is decomposed as g — ^r = i 9i^iA being the Pauli matrices . 

3 We take the d x 2_ y 2 order parameter as having its positive lobe along the the crystal a axis and its negative lobe along the 
&-axis. 1-1 
4 For symmetries obeyed by propagators and self energies, seeE3. 

5 It seems that eq. |^ must reflect itself in the boundary condition that should be used in a phgsomenological Ginzburg-Landau 
theory. It is unclear whether this has been taken into account properly in the literature (e.g.Lj) 
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i.e. the ratio is cutoff dependent. We must keep track of the the cutoff as soon as the T c -ratio is non-zero. The results 
given below are all at the cutoff n c (T/T c ) chosen equal to the integer part of 16T C /T. 

With the order parameter determined self-consistently, we can calculate the current-phase relation across a junction. 
For this we need the Green's function at the orifice. The propagator is computed along trajectories through the orifice 
(see Figure [j]). The boundary condition is boundedness far away on both sides of theJjiterface. At the orifice, the 
propagator is matched for continuity. An alternative route is the multiplication trickEd. This slight of hand takes 
advantage of exploding solutions along trajectories towards the junction. The matrix commutator of each pair of 
diverging solutions on the same trajectory at the pinhole (in fact an exploding solution and a decaying solution if 
viewed as propagating in the same direction) delivers the physical propagator in the orifice. The current density for 
a given phase difference x is calculated from the propagator in the orifice (at R = 0) by 

j s ( X ,T) = 2eN f v f Tj2 [ ^ cos 0Tr 2 [f 3 5 x (p, 0; e„)]. (8) 

Here Nf is the density of states at the Fermi surface, and > on the integral means that only the half sphere n • p > 
of directions are included. In general the Fermi velocity and the density of states depend on the position on the 
Fermi surface p. We assume a circular Fermi surface, for which the Fermi velocity is v/ = U/p. This assumption is 
not an unreasonable simplification. Buchholtz et a/.Ej compared calculations of the order parameter and the surface 
density of states using a cylindrical Fermi surface with the same calculations using a Fermi surface calculated from a 
tight-binding mode£3. They found very small differences between the two models. 

The phase dependent part of the Josephson energy density can be calculated from the current-phase relation as 

E a (x,T) = E s ( X o,T) + A f X d x 'j s {x',T). (9) 
Equation (^J) is valid for all the junctions discussed in the present article. 



III. THE TUNNEL JUNCTION 



As announced in the introduction, we wish to compare weak link junctions with tunnel junctions. A tunneling 
barrier is characterized by its transparency D(<p). Self consistent or not, the transparency of the junction may depend 
on the direction of the incident momentum (0 being the angle p makes with the surface normal). We have chosen a 
transmission coefficient proportional to exp(— 8 sin 2 (jj), i.e. peaked in the forward dir ectio n. This model transmission 
is used though out this paper with one exception. In the initial discussion in Section |V B| the transmission coefficient 
is chosen as all peaked along the junction normal. This we refer to as the forward (backward) tunneling limit. 



A. The self-consistent tunnel junction 



The self-consistent order-parameter calculation of Section II A works here as well since the tunneling potential 
barrier is high enough to perturt 



current we take the expression o: 



-the order parameters on the two sides only little. For calculating the tunneling 



UX,T) 



2eN f v f 



d(t> 

2^ 



D&) cos0 [g%(- X )g? (X) -$?0e)Si ("*)] 



(10) 



where gf"^ are the anomalous part of the propagator at the tunneling barrier (in our case at the surface) in either 
superconductor. 
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FIG. 2. Current-phase relation of a weak link with = 9n/40 and an = it/40 and temperature 0.1 T c (Bi g ). (a) 
self-consistent current relations for three different T c -ratios: (solid line), 10 -3 (dashed line) and 10" 1 (dotted line), (b) 
mimicking the self-consistent current-phase relation via a mixture of B\ g and i?2 9 . The choices are: Pure Bi g order pa- 
rameters in both superconductors (solid line), A R ' L (cj>) — As l3 — O.O5Ab 2b (0) (dashed line) and A L {<f>) = As l9 (</>) and 
A R ((j>) — AB lg (4>) — 0.05iAs 2s {(f)) (dotted line), (a) and (b) the dotted lines: states breaking time-reversal symmetry. 
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FIG. 3. Energy-phase relation for the same junction as in figure Q. Lines as in Fig. 



IV. NON SELF-CONSISTENT JUNCTIONS 



Most of the time we investigate the junctions with the order parameter self-c0nsisJ.en.tl3t calculated. As a point of 
comparison, we consider the case of a constant order parameter up to the surfaceEJ'tSEJ. This non self-consistent 
calculation can usually be done analytically. Comparing with the self consistent calculations we get a measure on the 
effects of surface pair breaking on the current-phase and the energy-phase relations. 

Consider a junction, weak link or tunnel junction, with order parameters left unaffected by scattering at the wall 
or the junction surface. Assuming tetragonal crystal symmetry, the bulk-gap function in d-wave superconductors is 

A Blg (p) = Acos2(0-a) (11) 

for order parameters of d x 2_ y 2 or B\ g symmetry and 

A B2g (p) = A sin 2(0 -a) (12) 

for order parameter of d xy or Bi g symmetry. A is the temperature dependent maximal amplitude of the gap and <j) 
the angle p makes with respect to the junction normal n. 

The expressions for the current through the junctions are calculated inserting the bulk propagators into equations 
and (|10|). The weak link— current is given byEil, and the tunnel-junction current is basically the well known 
Ambegaokar-Baratoff formulate. 
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A. Self-consistent pinhole versus non self-consistent pinhole 



We calculate and compare the current-phase (Fig. |J) and the energy-phase relations (Fig. ||) for a self consistent 
(a) and a non-self-consistent (b) weak link at different T c -ratios at T = 0.1T C . The crystal orientations are chosen 
as (Xl — 97r/40 and an = tt/AO. We see that even with a pure Bi g order parameter, the very similar current-phase 
relationships differ significantly from the sinusoidal form. The energy minima occur at some x, neither equal to 
nor 7r. If we admit, in a self-consistent fashion, a small £>2 S -order parameter, the critical current increases a little but 
the energy minima are not shifted in position. When, on the other hand, the T c -ratio is set equal to the temperature 
(dotted line in (a) and (b)), the system is just on the edge of the area in the parameter space where the superconducting 
state will spontaneously break time-reversal symmetry near the surface. There the current-phase relation no longer 
obeys the symmetry j s (x) = —jsi^x) an d the critical current is increased. 

Turning to the equally very similar energy-phase relations, we see that the symmetry of the time-reversal invariant 
relation has been lost at the edge of the instability against a time-reversal symmetry breaking state (TRSB). This 
applies also to the mimicked state in the non self-consistent case (dotted line in (b). The junction can in fact be in 
two distinct states (with different order parameters). The second state has E s (x) — ► E s (2n — x) (not shown). 

An overall conclusion is that the principal effect of pair braking on a surface reflects itself as a simple reduction of 
the critical current. 
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FIG. 4. Critical current vs. temperature in a d-wave/d-wave weak link at four different crystal orientations. Unit of current, 
2ttT c /R. 
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FIG. 5. Critical current vs. temperature in a d-wave/d-wave tunnel junction at four different crystal orientations. Unit of 
current, 2irT c /R. 



V. CRITICAL CURRENTS AND POSITIONS OF ENERGY MINIMA 



A. Critical currents 



In Figure |1| the critical current of a weak link between two d-wave superconductors at different crystal misorien- 
tations is shown as a function of the temperature. With the misorientation increasing towards ? the temperature 
dependence of the critical current evolves from linear in (1 — T/T c ) at an — to quadratic, (1 — T/T c ) 2 at cur = ? 
for the non-self-consistcnt (model) calculation a nd an even higher power for the self-consistent one. The solid line 
present in Figs.[|-|6| indicates again (see Section IV A) that the principal effect of the surface scattering is reducing 
the critical current. Adding a portion of a i?2 S component compensates for this reduction depending on the relative 
strengths of the two representations. At misorientations close to j, entering a TRSB state at low temperatures may 
boost the critical current. Below Tj RSB the temperature dependence of the critical current undergoes an abrupt 
deviation from the (1 — T/T c ) 2 behavior as clearly seen in Fig. [|(d). In other words, the temperature dependence of 
the critical current can serve as a detector of a TRSB-state at low temperatures. 

In Figure ^ the same temperature dependence as for the weak link in Fig. ^ is shown for a tunnel junction. Same 
as with the weak link, we see that the critical current is reduced by the surface pair breaking. Approaching olr = j 
extends the initial linear temperature dependence close to T c to lower temperatures. Right at an = ?, the current 
vanishes as long as the superconducting state is not a TRSB-state. 
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FIG. 6. Critical currents of a d-wave/d-wave junction (a-b) and of an s-wave/d-wave junction (c-d) vs. orientation angle az, 
at T = 0.1T C . Weak links, (a),(c); tunnel junctions, (b), (d) . 

The critical current at T — 0.1T C as a function of crystal orientation oil is shown in Fig. ^. Here an — 0, i.e, 
the a-axis of one of the superconductors is parallel to the junction normal. In Fig.|| a-b the junction is between two 
d-wave superconductors and in ^| c-d between a d-wave and an s-wave superconductor. Once more the TRSB-state 
has a most spectacular effect on the critical current. For the tunneling junction, the critical current no longer vanishes 
at a>L — § as in the time- reversal symmetric state. 

B. Positions of the energy minima 

We now turn to the positions of the energy minima. There is an essential difference between the angular dependences 
in weak-links and tunnel junctions. If we ignore the change in the order parameter near the surface, the forward 
(backward) tunneling limit of eq. (|l0|) is 

js(x) =iocos2ai cos 2or sin x. (13) 

This equation was hrst written down inS, and frequent use of it has been made in the literature (e.g.i'i). It can be 
rewritten as 

Js(x) = jo^{cos26 + cos2{2a L - 9j)sm X . (14) 

with 8 = oll — cur. Staring at the above (equivalent) equations, one observes that the current-phase relation is 
sinusoidal which leads to the energy minima of the junction always sitting at the values or 7r of the phase difference. 
The rest in the equations can at most change the sign of the current which cannot alter the conclusion that we are 
always dealing with a normal or a 7r-junction. As above the phase of the minimum of energy is referred to as Xmin- 
Fig. [?|a displays the areas in which Xmin takes each of its possible values in terms of the crystal angles. With pure 
Big this plot is correct at all temperatures also for the self-consistent solution. It is obvious, as well, that the critical 
current can only change sign where it vanishes, i.e., when oil or an equals 7r/4, 37r/4 etc. 
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a * 

FIG. 7. Comparison of angular ranges where tunnel junctions are 7r-junctions (in the forward-tunneling limit) between Bi g , 
(a); equation (^),(b). 



Now take a look at the pinhole, first at the Ginzburg-Landau limit, 
parameter up to the junction, equation (||) takes the form 

1 



At that limit and with a constant d. 



3b(x) = jo (cos 26 - — cos2(2a L -0))sinx 
15 



order 



(15) 



The argument above that leads to Xmin always equal to or it applies here as well. Because of the different coefficients 
in front of the second cosine terms, Eqs.(|l4|) and ( |l5| ) deliver very different zones of normal or ir behaviors, as 
highlighted in Fig. [?]. Equation (|l4|) assigns the junction orientation an overestimated weight in positioning the it 
behavior compared to equation (|15[) which emphasizes the misorientation of the two superconductors. This difference 
can be ascribed to the pinhole's allowing a larger incident angle for quasiparticles to be transported across the junction. 
Although illustrated with the non self-consistent order-parameter junction we shall make it clear later that the result 
is much more generak 

Chaudhari and LintJ measured the critical current from a hexagon of d x 2_ y 2 -superconductor to another crystal of the 
same superconductor in which the hexagon was inbedded. The misorientation between the two is oll ~ olr = f • The 
interfaces were the edges of the films which the authors interpreted as tunnel junctions. At the different boundaries of 
the hexagon the angles were equal to J,— -jg etc. (see Fig. |j) Eq. ( |l4|) gives the critical currents across the edges 
of the hexagon as j c = 0, ±0.433jo- This is in contradiction with the measured currents, which display a modulation 
around an average current. The psult of the experiment has been cited as evidence against the d x 2_ y 2 -pairing state 
for the cuprates (see for instanceEIJ) . Taking the crystal contacts for weak links, the results of Chaudhari and Lin do 
not rule out the d :c 2_ y 2-pairing state. In Fig. ^ the critical current at T = 0.1T C for various T c -ratios are plotted. The 
critical current is taken at the constant misorientation j . Picking the simplest case with no subdominant components 
the critical current oscillates around an average, just as in the experiment. The amplitude of the modulation is 
admittedly smaller than experimentally reported. Obviously the pinhole is a very idealized model for a weak link, 
and we cannot expect every small detail getting accounted for. 
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FIG. 9. Critical currents at a fixed misorientation ? but varying junction direction ql of a weak link at T — 0.1T C . The 
pattern has the period £ 



Away from the Ginzburg-Landau region, the positions of the energy minima were studied numerically. T = 
Q.lT c (Bi g ) was chosen. In Figure [lO] the position of the energy minima \m at a given crystal misorientation and a 
given junction orientation can be seen for the weak link. The corresponding graphs for the tunnel junction can be 
found in Fig. For comparison, the junction with a constant order parameter up to the wall without a subdominant 



component is included. Only one of the two phase differences (xin , Xm 



related through \ 



(2) 



2tt- 



- Xrn is displayed. 
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FIG. 10. Contour map of energy minima at T = 0.1T C in a weak link. Finite slopes are where the junc- 
tion is neither a 0-junction nor a 7r-junction. Full lines, self-consistently calculated results, dashed line in (a), con- 
stant-order-parameter-up-to-the-wall-model. pure _Bi 9 ,(a); T c -ratio 10~ 3 , (b); 0.1, (c); 0.2(d) . With the larger T c -ratios, 
the junction orientation acquires an increased importance as the amplitude of the Big -component depends on the interface 
orientation. 



The constant-order-parameter weak-link junction has zones where \m continuously evolves from zero to ir. A self- 
consistently determined order parameter leaves this picture unaltered, i.e. it does not lead to noticeable changes in 
the positions of the weak-link energy minima at zero T c -ratio. This is in contrast to the tunnel junction (not in the 
TRSB-state) where the changes continue being abrupt and the constant-order-parameter model and the self-consistent 
junction deliver very different pictures (see Fig. [ll] a) . In the weak link, Xmin is a strong function of the misorientation 
a L — can but depends only weakly on (or ocr) separately. The tunnel-junction Xm varies as a function of both 
the misorientation and the junction direction. This is true as long as the admixture of the i?2 S -component is kept 



small (in Fig. |10|.b T c -ratio is 10~ 3 ). Cranking up the T c -ratio (see Figs. lO.c-d.) adds structure to the boundary 
areas between normal and it behaviors of both junction types. Most structure is found when the superconducting 
state close to the junction breaks time-reversal symmetry with the mixture d x 2 +y 2 +id xy as the order parameter. The 
TRSB also smoothes out the jump between the zero junction and the 7r-junction areas in the tunneling limit. 

Higher temperatures only make narrower the zones of continuous change in \ m from zero to w. Above 0.5 T c this 
area has essentially vanished and the Ginzburg-Landau formula ( |l5| ) is an amazingly accurate description of the weak 
link. 
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FIG. 11. Contour map of the positions of the energy minima at T = 0.1T C for a tunnel junction. Full lines indicate 
self-consistent results, dashed line in (a) constant-order-parameter model, (a) pure Bi g superconductors (b) T c -ratio = 10~ 3 , 
(c) 0.1, (d) 0.2 (The difference between thffi^pH and dashed lines in (a) arises due to the existence of zero energy bound states 
near the surfaces in those orientations. seeE3t3) 



VI. EXPERIMENTAL IMPLICATIONS 



A. Small spontaneous fluxes at crystal interfaces 

An important implication of the present results bear on the interpretation of the experiment o!0. In this experiment 
Kirtley et al manufactured triangular co: hexagonal inclusions where the region inside was 45 deg misoriented with 
respect to the outside (similar to that ofl, see Fig. ||and Section 5.2). Fluxes spontaneously generated at the corners 
were measured with a SQUID. The fluxes found were usually small fractions of the flux quantum $ , neither integer 
nor half-integer multiples. Moreover, the crystal orientations themselves (with respect to the normal of the junction 
plane) seem to have no obvious correlation with the fluxes observed (Kirtley, private communication) . 

When the dimensions of the inclusion for the measurement of flux are much larger than the Josephson penetration 
length Aj, the flux generated at a corner should be %t{Xmin ~ Xmin) where Xmin are ^ ne position of the energy minima 
of the two interfaces meeting at that corner (see, e.g.E3). If one assumes junctions, such as the tunnel junction, which 
display Xmin = or ir as long as the order parameter does not break the-time reversal symmetry, then the experimental 
results force one to postulate TRSB at least at the grain boundarieafj. Even with TRSB, the observation that the 
crystal orientation themselves seem to have no significant effect on the fluxes observed would still find no explanation. 
If the bulk order parameter is pure B\ g then a TRSB order parameter near the interface can only occur for oil or an 
near 7r/4, whereas a.R — ul is not directly relevant (c.f. Fig. ^ a.) 

If the interfaces are weak-links, however, everything falls in place. Since Xmin is neither nor n in general even 
without TRSB (see Fig. |10|) . Our unit-transmission limit seems to be somewhat too extreme, however, in that Xmm 
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depends on — an but is almost independent on the individual ai or an. For pure Bi g , the theoretical fluxes are 
too small. This holds true also with a substantial Bi g component induced near the grain boundaries as long as we 
are not in a TRSB state (see Fig. In reality the grain boundary is more likely to have a smaller than unity 

transmission coefficient, so that Xmin would vary stronger with the a's (i.e, with pure i?i g , the portrait of Xmin would 
have to evolve from Fig. a to Fig. |^ a when the transmission decreases) . []. 

The condition oil — cur approximately equal to 7r/4|4s crucial. For example other experiments involving misoriented 
crystals,Bll3'tj (as well as the "frustrated" samples ofQ) all have misorientations close to (but sometimes not exactly) 
7r/6 and one of the a's is approximately zero, (or those related to these by symmetry) Our pinhole results for 
Xmin with the transmission coefficient equal to unity, differ only very slightly from the corresponding tunnel-junction 
values. ( For the junctions with qj, = 0, — ±7r/6, Xmin ~ 0.17-7T [| rather than x-m = for tunnel junctions without 
broken time reversal symmetry. ). The deviation for these new energy minima will be even smaller with a smallpj 
transmission coefficient and/or more directional tunneling. Thus our picture for the small spontaneously fluxes oft£l 
is not inconsistent with all these experiments. It would be interesting experimentally to try to observe the fluxes, if 
any, that nucleate at the corners of the inclusions with misorientation of ir/6 with judicious choice of a's. Obviously 
the results expected from the energy minima based on weak-links (Fig |Io| ) will be very different from those based on 
eq([l3|) (Fig |?] a) as well as those from tunnel junctions with TRSB near the interfaces (Fig |ll| c,d). 

B. The transmission coefficient of a junction 

J- pacpss-igrain boundaries have been measured. They range from ~ 10 3 A/cm 2 to ~ 10 5 A/cm 2 for T « T c ( 
e.g.Q,ffl,E3jij). A transmission coefficient can be inferred from these values. The critical current for a weak-link is 
dictated by the formation of phase-slip centers and is not equal to the thermodynamic current directly related to the 
transmission coefficient. RecentlytL 2 ] andlij have reported direct measurements of the Josephson penetration depth Aj 
via the observation of fluxes along the grain boundaries and the corner of the tricrystals (with a misorientation of 
~ 30 deg). We shall try to estimate the transmission probability \T\ 2 from their Aj. With A J 2 = ^?rw where w is the 
thickness of the magnetic field penetration layer along the barrier, and with the rough formula J c ~ irANfVfe\T\ 2 for 

the critical current density, as well as A^ 2 = AnN fv 2 e 2 / 'c for the London penetration length, we arrive at ^ ~ ^-|T| 2 

where £o = is the zero temperature coherence length. For YBCO, if we take Xl ~ 1400A , £o ~ 14A, with 
Aj ~ [im and w rj 2Xl, we get |T| 2 ~ 10~ 4 . Estimates for TIBaCuO are similar. At first sight then, |T| 2 seems 
small enough in the 3|0-deg misoriented films for one to assume the tunneling limit. If we accept the same order of 
mag nitude for \T\ 2 ir£9, then our scenario for the new Xm seem not very plausible. One needs to keep in mind, 
however, that the estimate is an average transmission coefficient. Electrical conduction may proceed through "hot 
spots" of much higher transparency. If we take the extreme limit of the pinholes with unit transmission, then the 
fraction of transparent area is « 1CP 4 . It is possible that the non-uniformity of the critical current densities cannot 
be seen in the scanning SQUID measurements if the distances between the pinholes are <C A,/. If d is the average 
pinhole-pinhole distance and (a/d) 2 ~ 10~ 4 , there remains reasonable-margin for the parameters a and d even with 
the restriction a << £0 << d « Xj for our calculation to apply. IntEl, the authors themselves emphasized that the 
fluxes would be much more localized than expected if Xj were ~ /xm, which may suggest that the boundaries are very 
non-uniform. 



VII. CONCLUSION 

In conclusion we have investigated the difference between high and low transmission junction in d-wave super- 
conductors, using pinhole as an extreme example in the high transmission limit. The effects of subdominant order 
parameter were also considered. 
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13 



(MF and SY), and Academy of Finland under contract No. 1081066 (MF) and on a sabbatical leave (JK). MF also 
acknowledges partial support from SFAAF, Abo Akademi and Magnus Ehrnrooths Stiftelse. Two of us (JK and SY) 
would like to thank Aspen Center for Physics, where this collaboration was initiated. 



1 Scalapino D.: Phys. Rep. 250, 329 (1995) 

2 Annett J., Goldenfeld N. and Leggett A. J.: preprint: [;ond-mat/9601060| , (1996) 

3 Wollman D. A. et al: Phys. Rev. Lett. 71, 2134 (1993) 

4 Brawner D. A. and Ott H. R.: Phys. Rev. B 50, 6530 (1994) 

5 Iguchi I. and Wen Z.: Phys. Rev. B 49, 12388 (1994) 

6 Wollman D. A. et al: Phys. Rev. Lett. 74, 797 (1995) 

7 Miller Jr. J. H., et al: Phys. Rev. Lett. 74, 2347 (1995) 

8 Chaudhari P. and Lin S-Y.: Phys. Rev. Lett. 72, 1084 (1994) 

9 Tsuei C. C, et al: Phys. Rev. Lett. 73, 593 (1994) 

10 Kirtley J. et al: Phys. Rev. B 51, 12057 (1995) 

11 Mathai A. et al : Phys. Rev. Lett 74, 4523 (1995) 

12 Tsuei C. C, et al: Science 271, 329 (1996) 

13 Kirtley J. et al: Phys. Rev. Lett. 76, 1336 (1996) 

14 Likharev K. K: Rev. Mod. Phys. 51, 101 (1979) 

15 Yip S.: Phys. Rev. B 52, 3087 (1995) 

16 Barash Yu. S. Galaktionov A. V. and Zaikin A. D.: Phys. Rev. B 52, 665 (1995) 

17 Ambegaokar V., deGennes P.G and Rainer D.: Phys. Rev. A 9, 2676 (1974) 

18 Buchholtz L., Palumbo M., Rainer D. and Sauls J. A.: J. Low Temp. Phys. 101, 1079 (1995); ibid 101, 1099 (1995) 

19 Mannhart J. et al , Phys. Rev. Lett 77, 2782 (1996) 

20 Kulik I. O. and Omel'yanchuk A. N.: Fiz. Nizk. Temp. 3, 945 (1977) 

21 Kopnin N. B.: JETP Lett. 43, 708 (1986) 

22 Kurkijarvi J.: Phys. Rev. B 38, 11184 (1988) 

23 Kurkijarvi J, Rainer D and Sauls J. A.: Can. J. Phys. 65, 1440 (1987) 

24 Fogelstrom M, Kurkijarvi J. and Sauls J. A.: Proceedings of the First International Conference on Quasiclassical Methods 
in Superconductivity and Superfluidity, Verditz, Austria. 

25 Matsumoto M., and Shiba H.: J. Phys. Soc. Japan 64, 1703 (1995); ibid 64, 3384 (1995); ibid 64, 4867 (1995) 

26 Zhang W.: Phys. Rev. B 52, 3772 (1995), Phys. Rev. B 52, 12538 (1995) 

27 Fogelstrom M., Rainer D. and Sauls J. A.: Phys. Rev. Lett. 79, 281 (1997); Fogelstrom M., Palumbo M., Buchholtz L. J., 
Rainer D. and Sauls J. A.: preprint 

28 Serene J. W., and Rainer D.: Physics Reports 101, 221 (1983) 

29 Palumbo M., Buchholtz L., Rainer D. and Sauls J. A.: Proceedings of the First International Conference on Quasiclassical 
Methods in Superconductivity and Superfluidity, Verditz, Austria. 

30 Sigrist M, and Rice T.M.: Rev. Mod. Phys. 67, 503 (1995) 

31 Thuneberg E. V., Kurkijarvi J. and Rainer D.: Phys. Rev. B 29, 3913 (1984) 

32 Radtke R. J., Ullah S., Levin K, and Norman M. R. Phys. Rev. B 46, 11975 (1992) 

33 Millis A. J., Rainer D. and Sauls J. A.: Phys. Rev. B. 38, 4504 (1988) 

34 Yip S.: J. Low Temp. Phys. 91 203 (1993) 

35 Ambegaokar V. and Baratoff A.: Phys. Rev. Lett. 10, 486 (1963). 

36 Sigrist M, and Rice T.M.: J. Phys. Soc. Japan 61, 4283 (1992) 

37 Tanaka Y. and Kashiwaya S., Phys. Rev. B53, 11957 (1996) 

38 Barash Yu S., Burkhardt H. and Rainer D., Phys. Rev. Lett. 77, 4070 (1996) 

39 Millis A. J.: Phys. Rev. B 49, 15408 (1994) 

40 Sigrist M, Bailey D. B. and Laughlin R. B.: Phys. Rev. Lett. 74, 3249 (1995) 



14 



